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Hamiltonian mixed systems with unbounded phase space are typically characterized by two asymptotic alge- 
braic laws: decay of recurrence time statistics (7) and superdiffusion (/?). We conjecture the universal exponents 
7 = /3 = 3/2 for trapping of trajectories to regular islands based on our analytical results for a wide class of 
area-preserving maps. For Hamiltonian mixed systems with bounded phase space the interval 3/2 < 7^ < 3 
was obtained, given that trapping takes place. A number of simulations and experiments by other authors give 
additional support to our claims. 
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Understanding the phenomenon of recurrence, introduced 
by Poincare in the end of the 19th century fP], has been a mat- 
ter of intense investigation over the last two decades in the 
theory of nonlinear dynamics. Let us consider the recurrence 
time probability P(r) defined as the probability of a trajectory 
to return at a time t > t to a pre-defined region. In the case 
of chaotic Hamiltonian systems the phenomenology of recur- 
rence is usually divided into two groups. For fully chaotic 
systems its long-time decay is exponential: P(r) ~ e~''^. 
This result can be understood recalling that: (i) the frequency 
of periodic orbits with period tq < t is given by the same 
relation, where h is the Kolmogorov-Sinai entropy [2J, and 
(ii) the Bowen's theorem guarantees the asymptotic equiv- 
alence of averaging over the phase space and over periodic 
orbits. The recurrence statistics of such systems is well un- 
derstood in the literature Nonetheless, most physically 
realistic systems are mixed, with the phase space consisting 
of chaotic and regular components. For these systems stick- 
ing to regular regions takes place, and the probability P{t) is 
believed to decay algebraically: 



P(r) 



(1) 



The algebraic decay ([T]) was first observed by Channon and 
Lebowitz jstl in the study of stochastic motion between two 
KAM surfaces in the Henon quadratic map. It is known 
that statistics of recurrence ([T]), clearly different from the fully 
chaotic one, is associated to anomalous transport 101, e.g., in 
a diffusive process where the mean square displacement in- 
creases with a power of t not equal to 1: 



(2) 



for an unbounded dynamical variable x{t). This fact has 
motivated the use of suitable statistical theories such as 
continuous-time random walk (CTRW) [8] and fractional ki- 
netics (FK) lyi] to describe such systems. An alternative ap- 
proach has been argued that the power law decay ([T]i is due 
to the hierarchical structure of phase space, which can be 
described by Markov tree models lid ll ill . Despite signifi- 
cant progress made with these theories, the understanding of 



algebraic decay ([T]) from first principles, e.g., from the mi- 
croscopic equations of dynamics, is still very limited. The 
purpose of this Letter is to understand the microscopic mech- 
anisms which determine the macroscopic algebraic laws ([T]) 
and (in through possible connections between them. 

For mixed systems the trapping of trajectories by the sta- 
ble islands gives rise to power law asymptotics not only for 
P(t), but also for the correlation function C(t). Karney [T^ 
showed, using simple arguments, that these two quantities can 
be related through C(r) c>c X^t^r ^(0- Therefore, if C(t) is 
asymptotically algebraic, then C(r) ^ tP{t) and P(r) sat- 
isfies Eq. ([T]i, since 7 > 1 guarantees the finiteness of mean 
recurrence times. Now, recalling the relation between the dif- 
fusion coefficient D{t), defined by 



(3) 



2r 

and the correlation function C(r): D{t) ^ J2t=i ^(^) 
If Eq. © is satisfied, then C (r) is also asymptotically alge- 
braic and we have the following relation 



7 + /3 = 3. 



(4) 



In short, by means of very simple and general arguments we 
show that, for sticking to regular regions in mixed systems 
with unbounded phase space (UPS), Eq. (HJ is true if, and 
only if, Eq. (|2]i is also true, since 7 > 1 and /3 satisfy the 
constitutive relation 

In order to find universal exponents for Hamiltonian dy- 
namics we will consider the following class of area preserving 
maps 

xt+i xt + Kf{Ot), 9t+i =0t+ xt+i mod27r, (5) 

defined on the cylinder —TT<9<Tr,~oo<x<oo, where 
K is the stochasticity parameter. Now, let us consider general 
UPS Hamiltonian mixed systems. Using the relation (|4|i, on 
basis of our analytical results for the system Q, and from 
numerical simulations and experiments by other authors we 
will conjecture that for such systems 



7 = = 3/2 



(6) 
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for the trapping of chaotic trajectories in the vicinity of is- 
lands of regular motion. We also confirm the Chirikov- 



Shepelyansky universal relation 111311 : 



7 = 3, 



13 = 0, 



(7) 



for sticking of trajectories in the vicinity of golden invariant 
torii, where (3 = corresponds to the null global diffusion for 
T ^ oo. Eq. ^ was conjectured in Ref. |13] on basis of 
numerical investigations of the paradigmatic Chirikov-Taylor 
standard map |6] at the critical parameter value K = Kc ~ 
0.97163540631. At this value the last invariant torus is criti- 
cal, e.g., susceptible to arbitrarily small perturbations. There 
is some controversy about the precise numerical observation 
of the 7 = 3 decay, but it is believed that this exponent should 
appear only for larger times [14]. If we reflect on this situation 
through Eq. (|4]i we can understand the difficulty of achieving 
a conclusive numerical answer: the KAM barrier at K = Kc 
is the frontier of two distinct transport regimes, namely fi = 
and /3 = 1. For the first case we have limT^oo D{t) = 
for K < Kc, while the second exhibits normal diffusion co- 
efficient D « 0.1(7^ - Kef for Kc < K < Kc + 0.1 L15J. 
To distinguish the two regimes on a scale of 10^ iterations for 
K — Kc ^ 10^^^ as done by Weiss et al. lll4ll it is necessary 
to detect ((A^a;)2) < 10"^^. 

To verify Eq. ^ we will first consider the standard map, for 
which J [9) = sin 6* and several studies of its transport prop- 
erties have been made f^i]. This map exhibits, in the vicinity 
of K = 21tt for I ^ integer, the following power law for the 
diffusion coefficient: 



D{t) ^ {\1\t) 



1/2 



(8) 



e.g., (3 = 3/2, obtained analitically in a recent Letter lfl6ll . 
Note that this power law is valid for all / ^ 0. We also observe 
that this result can be generalized for the wide class of maps 
(|5]l. If there is trapping of trajectories due to accelerator mode 
islands, which are created around a stable Q-periodic orbit 
{(a;t,6't)}i<t<D'then 



4+Q 



K 



Q 

E 

t=i 



f{0*t) = 2ln, 



(9) 



for all integer I ^ 0. The Q = 1 modes lead to enhanced 
diffusion, and the condition of stability for these is given by 
|2 + Kf'{0*)\ < 2. Let us now assume that there ai-e 9c 
such that f'{9c) = and f"{9c) 7^ 0. These conditions are 
required to ensure the existence of a relevant asymptotic form 
for Jo(e) = (2^)-i/:,rf^expK/(0)]: 

MO - r E + '^<^c/4]}, (10) 

C 

where Gc = (27r|/"(6'c)|)"^/^ and <j)c = sgn[/"(6lc)]. Con- 
sider now 9* = 9c + e. We can always choose e small enough 
such that -4:/\Kf"{9c)\ < sgn{K)(f>c€ < 0, in accordance 
with stability condition, and Kf{9c) 21tt as in Eq. (|9]l 



for (5 = 1, since f{6c) ^ BlTll . Fortunately, by the same 
arguments showed in Ref. fl^, we also have 



DiT)^Kj2MjK) 



\r) 



1/2 



(11) 



otherwise the sum in Eq. ( fTTI ) gives terms of the type 
J2li=i exp(zg77) / whose sum converges by the Cauchy in- 
tegral test for T ^ 00. The effective contribution of Q > 2 
modes for enhanced diffusion is subject to several constraints. 
In addition to meeting Eq. (|9]l, the stability condition of such 
modes follows from the trace of the monodromy matrix 



n 



Kf'{9;) 



1 i + Kf'i9;) 



(12) 



All of these restrictions make the modes less significant as Q 
increases. In short, the algebraic law (fTTI) is closely associated 
with the existence of critical points 9c for which the asymp- 
totic form (fTOb exists. We point out a strong fact supporting 
the universality of f3 = 3/2 (and thus 7 = 3/2): this exponent 
does not depend on the general form of f{9) nor the number 
of critical points, just the existence of a single critical value 9c 
is enough to justify it. 

In order to reinforce the argument of universaUty for Eq. ^ 
we can mention a number of simulations and experiments that 
also point towards the universality of these exponents. Karney 
1 12|] computed P{t) for the quadratic map and Chirikov and 
Shepelyansky jisll did the same for the whisker map. Their 
numerical results agree with 7 = 3/2. In a recent Letter, 
Cristadoro and Ketzmerick 11 ill conjecture the universality of 
only Eq. ([TJ. They used a Markov tree model, originally de- 
veloped by Hanson et al. [10], with random scaling factors for 
the transition probabilities. The long-time algebraic decay ([T]) 
was then obtained. In order to check the universality of the 
Eq. ([1]) they studied numerically an specific area preserving 
map of the type (|5]l (for which 9c exists!) for several combi- 
nations of parameters. The averaged exponent obtained was 
also very close to 3/2, namely 7 = 1.57. There are also in 
the literature a number of experiments for the study of 7 and 
/3 from the perspective of CTRW approach, for which Eq. (HI 
is also valid. Zumofen and Klafter [8] studied the Levy-walk 
statistics of standard map for two specific values of stochas- 
ticity parameter: K = 6.716809 and K = 6.476939, related 
to Q = 3 and Q = 5 accelerated modes, respectively. The 
exponents obtained by them are: 7 = 1.8 and (3 = 1.2 for 
Q = 3, and 7 = 1.6 and /3 = 1.4 for Q = 5. Note that 
K = 6.476939 is closer to K = 2tt, for which the power law 
diffusion is more precise iflill . This explains the best fit of the 
Q = h mode with the universality assumption (|6]l. Sanders 
and Larralde [19] studied the ocurrence of anomalous diffu- 
sion in polygonal billiard channels. They found 7 = 1.58 and 
/3 = 1.81 for a particular billiard model. A similar experi- 
ment was also performed by Schmiedeberg and Stark [20] for 
an extended billiard model with honeycomb geometry, get- 
ting (3 = 1.72, and by Dettmann and Cohen [21] for a wind- 
tree billiard model obtaining 7 = 1.4. We should observe 
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TABLE I: Algebraic exponents for some UPS Hamiltonian systems. 



SIMULATIONS - EXPERIMENTS 


7 


/3 


quadratic map [6. 12] 


1.44 


- 


whisker map [6, 18] 


1.44 


- 


standard map for K = 6.476939 [8] 


1.6 


1.4 


parainetric map model 111] 


1.57 


- 


continuous time Hamiltonian model ]22] 


1.55 


1.4 


flow with hexagonal symmetry 17] 


- 


1.4 


two-dimensional rotating fluid flow ]23] 


1.6 


1.65 


fluid flow with no KAM surfaces ]24] 




1.55 


zigzag billiard ]19] 


1.58 


1.81 


polygonal Lorentz billiard 119] 




1.4 


honeycomb billiard ]20] 




1.72 


wind-tree billiard ]21] 




1.4 


mean values 


1.54 ± 0.07 


1.53 ± 0.16 



that, usually, there is a coexistence between normal diffusive 
(3=1 and superdiffusive 1 < /? < 2 processes. Thus, the 
algebraic laws ([1]) and (|2]i cannot be accounted for quanti- 
tatively by purely CTRW approach unless sufficiently large 
times are reached. In this sense, Gluck et al. [22] studied the 
process of chaotic diffusion for the continuous time periodic 
Hamiltonian H^{p^q,t) = + cos^(q) cos^((jjt), where 
— oo < g < oo. In order to include the normal diffusive ef- 
fects they used a convolution of normal and Levy propagators 
summed over all Levy flight times. For cj = 0.8 they found 
7 = 1.55 and (3 ~ 1.4, in good agreement with Eq. Fi- 
nally, a numerical study of a stationary flow with hexagonal 
symmetry |j7t] and two experimental studies on chaotic trans- 
port in fluid flows lEsi |24 ] also agree with 7«/3«3/2orat 
least with /3 « 3/2. All results described here are summarized 
in Table U As we can see, the mean values obtained for 7 and 
(3 are in agreement with the universality of Eq. (|6]l. 

An interesting connection between the Eq. ^ and the hi- 
erarchical structure of p hase space can be found in a work 
of Zumofen et al. [25']. They study the motion of a parti- 
cle in layered medias through its random velocity field along 
the X direction accompanied by normal diffusional motion in 
the space y transverse to x: x{t) = v[y{t)]dt. One of 
the transverse spaces studied by them is a system with hierar- 
chically arranged connections between layers, very similar to 
Markov tree models used in Refs. iiSinil- Assuming a power 
law decay for the correlation of velocities Cu(t) ^ r^^ 
they conclude that (3 = 2 — min(A, 1). By Eqs. ([T]l and 
C(r) ^ tP{t) we have A = 7— 1, which confirms once again 
7 + /? = 3, since 1 < 7, /? < 2. They also studied another 
model of walks whose transverse space is one-dimensional, 
finding /3 = 3/2. However, the statistical assumptions used in 
this model do not correspond to the dynamics of the map (|5]l 

So far, our results include only UPS Hamiltonian mixed 
systems, for which global diffusion takes place. Nonethe- 
less, an additional observation on Hamiltonian mixed sys- 
tems with bounded phase spaces (BPS) is in order. Suppose 



that for times r > Tq the algebraic asymptotic regimes hold: 
Pu{t) = CuT~^^ and Pbir) = CbT^^*, for unbounded (u) 
and bounded (b) systems, respectively. On the one hand, ac- 
cording to Poincare's recurrence theorem [1], the return of al- 
most all trajectories in a finite time is assured only for BPS 
systems. On the other hand, 7„ = 3/2 > 1 guarantees the 
finiteness of mean recurrence times. Therefore, we also have 
finite return times for almost all trajectories in UPS systems. 
Let us focus now our attention on the trapping of trajectories 
by regular islands for both systems. As it is typical in such 
scenarios, the phase space is composed by islands of stability 
immersed in a sea of chaotic trajectories with no KAM sur- 
faces. Occasionally, there may also be nonhierarchical bor- 
ders between regular and chaotic regions for some types of 
BPS systems known to be sharply divided phase space ll27ll . 
Now, let us construct for each arbitrary BPS system with 
phase space i7f, its UPS periodic extension defined through 
translations x{t) x{t) + £. Here, ^ is a vector of the Bra- 
vais lattice for which fti, is the first Brillouin zone and x{t) 
becomes unrestricted. Such extension makes long recurrence 
times more likely. Thus, comparing the process of recurrence 
in both cases we must have Pu{t) > Pb{T) which implies, 
through universality of Eq. (|6]l, that 7b > 7u = 3/2 for 
aU T > max[ro, (Cb/C„)i/7.-3/2]^ = = 3/2 if 

Cb < C«- For sticking of trajectories in the vicinity of golden 
invariant torii a similar comparison is valid. In both cases the 
trajectories are equally bounded by KAM surfaces, so that we 
have 76 = 7u = 3 by Eq. (|7]l. Thus, the maximal trapping for 
UPS systems should correspond to an upper limit for the max- 
imal trapping in BPS systems. Finally, we have the inequality 



3/2 < 7b < 3. 



(13) 



The result dlSll has been confirmed by numerous computer 
experiments ITA 12811 . Special attention should be be payed to 
Hamiltonian systems with sharply divided phase space, where 
KAM torii coexist withergodic components of positive mea- 
sure. Altmann et al. li27ll studied numerically these systems 
(mushroom billiards and piecewise-linear maps) and their re- 
sults suggest the universality of the exponent 7b = 2 in the 
phase space ergodic components since hierarchical borders 
are not present. 

In conclusion, we show that the distribution of recurrence 
times and the power law diffusion for UPS Hamiltonian mixed 
systems are two macroscopic algebraic laws intimately con- 
nected through their asymptotic exponents: 7 + /3 = 3. This 
relation confirms the Chirikov and Shepelyansky conjecture 
of universality 7 = 3 for sticking of trajectories in the vicinity 
of golden invariant torii. We also conjecture (our main result) 
the universality of exponents 7 = /3 = 3/2 for trapping of 
trajectories in the vicinity of stable islands. This result was 
based on the analytical calculation of (3 = 3/2 for a wide 
class of area preserving maps. It therefore represents a sig- 
nificant step towards elucidating the mechanism of the direct 
extraction of characteristic exponents in generic Hamiltonian 
systems. For trapping of trajectories to regular regions in BPS 
Hamiltonian mixed systems, the interval 3/2 < 7b < 3 was 
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established, given that trapping takes place. Finally, our re- 
sults agree with numerous simulations and experimental re- 
sults existing in the literature. We should remark, however, 
that because the algebraic regime in Hamiltonian systems can 
take an arbitrarily long time, in general, it is not even clear 
for some systems whether the 7 = = 3/2 was sufficiently 
reached. Another important aspect concerns the usual lack 
of exact power law formulas for superdiffusion coefficients: 
the non-acquaintance of their parametric dependence may sig- 
nificantly affect numerical estimates of such exponents. As 
previously mentioned, the 7 and (3 values obtained in Ref. 
(SD for the standard map with K = 6.476939 are closer to 
the foreseen values 7 = /3 = 3/2 than those obtained for 
K = 6.716809. This occurs because the algebraic law Q ap- 
plies only to values of K very close to 21tt, since f{6c) = 1 
for such case. Now, let us return to the Gluck et al simulations 
1 2211 . They conclude that f3 « 1.4 obtained for uj = 0.8 is not 
universal because of the value P « 1.33 found for uj = 0.88. 
However, the Hamiltonian H^{p, q, t) studied by them can be 
derived from Hk{p, q, t) — + K cos'^{q) cos^(t), since 
the following scaling transformations are made: H H/ A, 
p — > p/V\, q ^ q, t \/At, and K ^ K. Choosing 
A = we have K — uj^'^l'^ relating Hk and H^. The 
Hamiltonian Hk has a similar structure to the kicked rotor 
Hamiltonians that lead to map (|5]l lEoll . Therefore, as dis- 
cussed above on the standard map, it is not surprising to note 
some discrepancy between the two values of ji blindly fitted 
for the two u: values, noting that = K^^I^Dk for the 
momentum variable p. 

The author thanks A. Saa, W.F. Wreszinski, C.J. A. Pires, 
and R. Vicente for helpful discussions. 
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